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CONVECTION ELECTRIC FIELDS AND THE 
DIFFUSION OF TXUPPED 
MAGNETOSI>HERIC RADIATION 
Thomas J. Birmingham 

ABSTRACT 

We explore here the possible importance of time-dependent convection 
electric fields as an agent for diffusing trapped magnetospheric radiation in- 
ward toward the earth. Using a formalism (Birmingham, Northrop, and 
Falthammar, 1967) based on first principles, and adopting a simple model for 
the magnetosphere and its electric field, we succeed in deriving a one dimen- 
sional diffusion equation to describe statistically the loss free motion of mir- 
roring particles with arbitrary but conserved values of the first two adiabatic 
invariants M and J. Solution of this equation bears out the fact that reasonable 
electric field strengths, correlated in time for no longer than the azimuthal drift 
period of an average particle, move particles toward the earth at a rate at least 

an order of magnitude faster than electric fields whose source is a fluctuating 

/ 

/ 

/ 

/ 
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current on the magnetopause. 


CONVECTION ELECTIUC I^IELDS AND THE 


DIFFUSION OF TRAI^PED 
MAGNETOSPHEKIC RADIATION 


INTRODUCTION 

There has been growing rapidly a body of literatta'e concerned with the 
theory of interaction of magneto spheric particles with collective oscillations of 
the magneto spheric medium (Wentzel, 1961; Parker, 1961; Dragt, 1961; Dimgey, 
1963; Cornwall, 1964; Chang and Pearlstein, 1965; Kennel and Petschek, 1966; 

Cocke and Cornwall, 1967; Roberts and Schulz, 1968). The thesis of this work 
is that at heights where interparticle collisions are infrequent, wave-particle 
interactions are an effective agent for depleting the energetic trapped population: 
successive wave induced scatters of a particle’s pitch angle cause it to mirror 
at altitudes increasingly closer to the earth, until ultimately the particle reaches 
the collision dominated level of the ionosphere where its mirror motion is 
disrupted. Impressive corroboration of esqDeriment with theory is presented by 
Kennel and Petschek (1966) who show that for a range of L-values (>4) equatorial 
fluxes of both energetic protons and electrons lie near values expected theoretically 
when particles and waves are in qua si-linear equilibrium. 

The existence of steady-state particle fluxes in the presence of such losses 
obviously requires a balancing input. We here consider in detail a version of the 
theory that particles are injected at tlie magnetopause and then diffuse in L, the 
equatorial crossing distance of a magnetic field line expressed in numbers of 
earth radii R^. 


1 





Diffusion in L is likewise a topic which has received extensive attention. 

An exhaustively explored (Kellogg, 1959; Parker, 1960; Davis and Chang, 1962; 
Tverskoy, 1964; Dungey, 1965; Nalmda and Mead, 1965; Conrath, 1.967) idea has 
been that L-diffusioa is driven by a time varying electric field which arises 

IT 

from fluctuations in the currents separating the terrestrial and interplanetary 
magnetic fields. The calculated particle diffusion rate is generally sufficient 
to bJilance interparticle scattering losses (Nakada and Mead, 1965; Newkirk 
and Walt, 1968) but appears grossly inadequate when wave particle interaction 

losSjCS are taken into account (Kavanagh, 1968). 

i 

iCladis (1966) has proposed a slightly different mechanism: at very low 

i 

L“values ionospheric current variations produce fluctuating electric fields which 
drive the diffusion process. The effectiveness of such a mechanism diminishes, 
however, as a function of altitude in the magnetosphere® 

Somewhat in contrast, one of Cornwall’s (1968) speculations is that L-diffusion 
is a magneto spheric manifestation of Bohm diffusion (1949), an anomalously large 
diffusion rate across magnetic field lines observed in many laboratory plasma 
ekperiments (e.g., cf. Pe?se et aL, 1966). It is thought that this enhancement 
above the level of classical, colli si onal diffusion is the result of particles inter- 
acting with turbulence in the plasma. The turbulence in turn is a consequence of 
the fact that laboratory plasmas have density and/or temperature gradients and 
hence are susceptible to the growth of unstable drift waves (cf. Kadomtsev, 1965). 
In the magnetosphere similar gradients exist, and it is possible that Bohm dif- 
fusion is operative, particularly in the vicinity of the plasmapaiise. 
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Oul’ own viewpoint is in line with the thinlcing of Falthammai% who has con- 
tinually stressed the importance of electrostatic fields (V.>, E = 0), varying in 
time and of magneto spheric spatial dimensions, in the dynamics of trapped 
radiation. Specifically, we consider as the driving force in L-diffusion the 
electrostatic field 

Vf >; B 

E(ri t) “ VV **■ (1) 

associated with the ideal hydromagnetic flow of low energy magneto spheric 
plasma. (Oui* demarcation between low and high energy particles is roughly 
5 keV, an energy at which gradient and curvature drifts become important com- 
ponents of the guiding center motion at magneto spheric heights for tj^ical 
electric fields. Low energy particles essentially follow the hydromagnetic con- 
vection pattern and are of interest to us only as the source of the polarization 
electric field.) Falthammar (1965), Brice (private communication), Obayashi 
and Nishida (1968), Kavanagh (1968), and Cornwall (1968) have each suggested 
the possible importance of the convection electrostatic field to the diffusion of 
energetic particles. 

In this paper we adapt the diffusion theory developed from first principles 
by Birmingham, Northrop, and Falthammar (1967) to a simple model of the mag- 
netosphere and its impressed convection electric field. No restriction to equa- 
torial particles is made. The one dimensional diffusion equation which statis- 
tically describes loss-free particle motion in this model is then solved. Results 
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are compax’ed with those frorn. a model in which magnetopause currents are the 
source of the magnetosphcric electric field. We conclude tliat convection electric 
fields move particles radially at a rate an order of magnitude faster than the 
electric fields arising from magnetopause currents. 

Unless otherwise specified, Gaussian (cgs.) units are used throughout this 
paper. 

MODEL 

We consider the following elementary model of the earth’s magnetosphere. 
The magnetic field B, constant in time, is due solely to a magnetic dipole of 
moment /xat the origin of an r, t?, ^ spherical coordinate system. In the con- 
ventional sense, m is aligned with the south normal to the ecliptic plane, co- 
latitude t? is measured from the axis, and azimuth <t> is measured clockwise 
from the antisolar meridian (i.e., 0 " tt/ 2 on the dawn meridian). Modification 
of the magnetic field by either the energetic particles or the cold plasma com- 
ponent is neglected in this low-/S model. 

The electric field E is irrotational, variable in time, and (from Equation 1) 
everywhere perpendicular to B. It has the fea,ture, essential for L -diffusion, that 
it is asymmetric about an energetic particle’s longitudinal drift path in the dipole 
field. We describe E by the potential V 


-«r 




n 

ft 


A being a positive, time -dependent amplitude. The form Equation (2) is the 
fundamental (m - 1) asymmetric mode in FMiammar’s (1965) Fourier expansion 
of a general longitudinally dependent potential. Since r sin"^ d and <p are both 
constant on dipole field lines, B lines are equipotentials and E • B is zero. In 
the d “ 7//2, equatorial plane of the dipole, the electric field derived from Equa- 
tion 2 has the magnitude A and is uniformly directed, dawn to dusk. 

In this simple model, low energy plasma in the equatorial plane flows in 
straight lines (oriented along " tt) from the night to the day side of the mag- 
netosphere. The model thus crudely represents (in a mathematically tractable 
manner) the flow pattern as depicted by Levy et al. (1964) for Dungey's (1961) 
field line merging model of the magnetosphere (cf. also Petschek, 1964; Brice, 
1967). Except for regions near the magnetopause our model flow is also crudely 
representative of the pattern envisaged by Axford and Hines (1961) in their closed 
magnetosphere model. (We argue further that by allowing A in Equation 2 to 
reverse sign beyond a certain large L = Lj^, an even better approximation to the 
Axford and Hines picture is obtained. The diffusion of energetic particles in our 
analysis depends only on the autocorrelation of A. Thus, as long as the statistical 
properties of A are the same for all L, our results are valid even for L > Lj^.) 

We treat A as a fluctuating quantity with an average magnitude of 4 x 10“'^ V/m, 
a value typical of hydromagnetic models of the magnetosphere. The fluctuations 
reflect similar variations in the intensity of the solar wind, to which the mag- 
neto spheric flow pattern is directly coupled. A one hour time scale seems 
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(Oarponter aiicl Stone, 1067; ObayasM and Nishida, 1068; Bi’ice, private com- 
munication) appropriate for fluctuations in A(t). Wo shall find that the mean 
square fluctuation amplitude G, - and the autocorrelation time of the 

fluctuations are quantities which scale the diffusion time of the energetic particles. 

Our simple model takes no account of the electric fie- 4 associated with the 
rotation of a finite radius, conducting earth. In hydromagnetic models of the 
magnetosphere, the rotation electric field is considered to dominate the dynamics 
of low energy plasma at radial distances as far out as the plasmapause. Some 
evidence exists (Nishida, 1966; Carpenter and Stone, 1967), however, that the 
pattern of outer magnetosphere electric fields persists, though perhaps in a 
recessive. role, in the region dominated by rotation. The electric field associated 
with the rotation of the earth is both nearly longitudinally symmetric and nearly 
constant in time. (The tilt of the earth’s magnetic axis with respect to its spin 
axis introduces a 24 hour periodicity into the rotation electric field. Such a time 
variation is too slow to be consequential in the diffusion of energetic particles.) 

For L -diffusion the rotation field plays no significant role. 

We realize fully that the model of electric and magnetic fields adopted in 
this paper lacks the detail of extant qualitative hydromagnetic models. However, 
for our purposes — • an analytical study of L-diffusion •— we feel that the present 
model exl;racts the salience of its more complicated counterparts and simultan- 
eously affords a simplicity necessary for achieving mathematical results. 
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THE DIFFUSION EQUATION FOB ENEBGETIC PAKTICLES 

For onorgetic particles in our model, the ratio of particle gyro-radius 
to the scale length € (tyt)ically IRp ) of electric and magnetic fields is small; 

t f 4L3 /mT y 10^ 1 . 


Here W is the energy (in eV) of an energetic particle of mass m (in grams) at a 
distance LRp from the earth. Further, the typical one hour time variation of the 
electric field lies on the slow time scale of the azimuthal drift of these particles 
in the dipole magnetic field. Finally we note that the ratio E/B of the field mag- 
nitudes, electric to magnetic, is at least 0(e) in smallness. 

Northrop’ s (1963) systematic development of tlie adiabatic motion of charged 
particles asserts that for the conditions stated in the preceding paragraph, 
energetic particles move in such a fashion that the first two adiabatic envariants 
of their motion — the magnetic moment M = mvxV^2B and the longitudinal in- 
variant J " m ^ V|| ds — are conserved. (Here i and ll deact*. velocity components 
perpendicular and parallel to the dipole magnetic field, and the s -integral is 
extended over a complete guiding center bounce path.) The third or flux invariant 
«l> is not preserved in the present treatment. 

Magnetic moment conservation implies that particles moving inward in the 
dipole magnetic field gain Idnetic energy. A rough estimate based on both M and 
J conservation 


pyjl/2 

J 


const a 


(''lOoq L 


L ^/2 ctnXgq 


(4) 
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also indicates that inward moving particles flatten their equatorial pitch angles 
K , In doing so, such particles contribute to a pitch -angle anisotropy which can 
drive microinstabilities and lead to diffusion into the loss cone of the magnetic 
mirror, 

All particles with the same values of the invariants M and J and with guiding 
centers lying on the same dipole field line (we here identify a dipole line by the 
constant values of the Euler potentials, /3 = (^, a = - sin^ d/r “ - /VLR^, on it) 
have practically the same experience over the course of one bounce, since the 
electric field is changing on the much longer drift time scale. In our treatment 
we consider a number of (identical) boimce-sveraged guiding centers with invar- 
iants M and J as closely representing the actual M, J particles gyrating about and 
bouncing along the line a,/3. Equations for such bounce averaged guiding centers 
were first obtained by Northrop and Teller (1960), who performed a bounce 
average of the ordinary guiding center equations of motion. 

Within the framework of bounce-averaged guiding center theory our treat- 
ment of the energetic particles is a statistical one: we adapt to the present 
model a general theory (Birmingham, Northrop, and Falthammar, 1967, here- 
after referred to as BNF) for the diffusion of guiding centers, the diffusion in 
the BNF theory being due to small electromagnetic fluctuations which conserve 
M and J. Only statistical propertie.s of the fluctuations are assumed to be known. 

In the present model, of course, the magnetic field is static and the BNF treatment 
simplifies greatly. 


8 



I 

H 


III Appendix A adaptation of the BNF theory- to the present model is carried 
out and the equation 


r?<Q>(a. MJ,t) 



(5) 


is derived. Here <Q)’ is the /3-average of the guiding center density <Q)>(a, /3, M , 


J, t) defined for a four dimensional ~ a, /3, M,J — phase space: 


f>2'rr 

<Q> = <Q> 

Jo (6) 


itself is a statistical average, this average (denoted by < ]» being carried 
out for the same set of M,J particles over an ensemble of fluctuation systems. 
The /3-average is carried out at constant a, M , J, and t. Use of the definition 
a = - M/LRg and the relationship 


dL ^ /7\ 

= aj n(L.MJ.t) = 

leads to the following more familiar form of Equation 5, 

Jt ~ A, > (8) 

where 

Re^ 

f(L) = (L) . (9) 

r' 
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Wliile Equation 8 is identical in form with the equation generally introduced in 
L~diffusion studies, it is here valid without the restriction to J = 0 particles 
inherent in most past work. 

Were there no electric field, M, J particles would drift longitudinally at 
constant a in the dipole field with an angular drift frequency o)jj (a, M, J), Be»- 
cause of the fluctuating electric field, however, suchM, J particles drift across 
a, changing tiieir kinetic energy in the process. Equations 5 and 8 describe this 
diffusion across dipole field drift shells, all gyro-phase, bounce-phase, and 
longitude dependence having been averaged from the problem at the present level 
of descrij)tion. 

Northrop (private communication) has shown the following relationship be- 
tween (a, yS, M, J, t) and the particle Eux j differential in kinetic energy W 
and solid angle H; 

<Q> ~ ■2'W ' (10) 

Since is independent of arc length s along each field line, j/W is similarly 
s-independent. At each value of s along a line a, /3, however, the flux j is to be 
measured for a different solid angle H and a different energy W, n and W being 
determined by the conservation of M and J on that field line. 


10 


The diffusion coefficient is derived in tei'ms of SV, the fluctuating part 
of the potential, in a straightforward manner from BNF results: 



Subscripts here indicate the functional dependence of the quantities involved; thus 
is a fimction F(a, t) while is F(af /5 - r, t- 

Using Equation 2 and the definition of a and /S we further obtain 






dr § A( t - 7 ) oA( t ) cos 




( 12 ) 


Dependence of on M and J for this case where magnetic field lines are equi- 
potentials occurs only through co^ 


J) " 


2c 

— MB^ (a, M, J) 


J 


2 4MB T(a,M, J) 


(13) 


Equation 13, a form due to Northrop (1966), gives for a particle with invariants 
M and J whose guiding center is on dipole line a at any /3, in terms of M, J, a, B^, 
and T. Here and T are respectively the mirror magnetic field strength and 
the bounce period of such a particle. 
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In order that the longitudinal invariant J exist, E must be at least 0(fe) in 
smallness compared witli the 0(1) magnetic field; as in BNF, we consider here 
the yet more restricted case where E is 0(e^ ), e « S « 1, With this ordering 
and the form of the coefficient as given by Equation 12, our diffusion equation. 
Equation 5, is correct through 0(eS^]. 

In passing we note that 



so that, using Equation 12, 


/ (AL)A 

\ / 


L’Re'* 



(14) 






*•00 

dr <8A(t-r) SA(t)> 

0 


cos co^ r . 


(15) 


Equation 15 exhibits a familiar dependence (F^thammar, 1965; Cornwall, 1968) 
on a) the sixth power of L (explicitly) and b) the value of the power spectrum of 
the fluctuating part of tb.e equatorial electric field at the resonant frequency 

0) COj^m 

A reasonable direction to proceed, in view of the paucity of direct experi- 
mental evidence of electric fields and their time variations, is to assiune that the 
autocorrelation <(8A( t - r) SA( t))> has the form 


<SA(t -r) SA(t)> 


d exp - 



(16) 
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stationary in time (independent of t) » This choice is similar to one made by 
Cornwall (1968) and is here motivated by the idea that the equatorial electric 
field A exhibits no periodic amplitude variations, is always directed from dawn 
to dusk, and is random on the time scale on which the solar wind executes time 
variations of large spatial extent. (The correlation time r is thus typically one 


hour.) 


Substituting Equation 16 into Equation 12 and integrating, we obtain 


D 


aa 


4a 2 


fn T (1 exp 




(17) 


Given a value of r^, note that for fast particles with short drift periods 
Tp = 277/^1 u>j^ I « is e^onentially small . Such particles experience an es- 
sentially time independent ensemble average electric field which drives them 
to smaller a for one half of their azimuthal drift period and to larger a for the 
other half; over a complete drift period the two effects cancel. 

For slower particles, , we can represent approximately as 


D 


aa 




yiT r a 


(18) 


independent of a particle’s M and J . In the following section we shall consider 
the solution to the diffusion Equation 5 with Equation 18 as the diffusion coef- 
ficient. Note from Equation 13 that the approximation = 2rr/\oi^ \ » restricts 
M , J particles which can be effectively diffused to an energy range determined 
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If* 


n 

if 


ff 


from 


7T \q\ II 1 

= "w n 

wii w/ 


» T 


c 


For given W andL, the greatest restriction is on J = 0 particles. 


(19) 


SOLUTION OF THE LOSS-FKEE DIFFUSION EQUATION 

Pursuant to the discussion in the previous section we now write down the 
equation 

a<Q> _ a 

dt " “c^a 

which describes the loss free diffusion of particles whose conserved M and J 
satisfy (19) at all a. 

The general solution to Equation 20 regular at a = 0 (L = <») and a = - oo(L = 0) 
for the initial condition t = 0) = 'K<^) can be worked out in a straight- 

forward manner. We refer the reader interested in the mathematical steps 
leading to this solution to Parker’s (1960) paper and here merely quote the result 




L\ 4a2 


r d 

■ C j 


^<Q> 

da 


( 20 ) 




dy y3/2 J 


3/4 


xy^ 

"T* 




( 21 ) 
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In Equation 21 J3/4 is the Bessel Fimction of order 3/4 and r " 'fn dt/4 

is normalized time. Note tliat r depends linearly on both the correlation time 
and (1, seen from Equation 16 to be the (time independent) mean square 
fluctuation ^[sA(t)] ^ of A. 

The choice of 0(a) " N5 (a “ a^), corresponding to a initially spiked at 
L. ~ “/i/a. Rg, facEitates integration of Equation 21. The time development of 
Q)>is then essentially tlie same as if 0were a Gaussian whose dispersion <r^ is 
much smaller than L.^. 

Using the S-function form of 0, we pei’form the integrations in Equation 21 
and obtain 


<Q> = N 




a’ + a’ 

^ I 


16r / -^3/4 V 8r 


a.^ 

1 


13^4, the 3/4 order Bessel function of imaginary argiunent, has the simple integrt 
representation (Gradshteyn and Ryzhik, 1965) 


^3/4 ~ 


iVrfi] I. ' 


l_jj2^i/4exp (zx) dx , 


r being the usual gamma function. With 13/^ expressed in integral form and 
identified as /VtRg, <C Q/has the form 


<Q> 


32r(| r(4) r’/-* 


(J- + ' 

16Rg'^r\L.'^ Vj 


dx (1 - exp 


8 L 2 l 2 r 
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The L, r dependence of has been investigated numerically for a case in 
which the initial position of the 8-function is L. = 8, The results appear in 
Figure 1: is in arbitrary units and T ” 2.4 L.'' r = l,4xl0®(r^(l)t 

[t and are inhrs., (X in Values fx/R^ = .3 gauss, L. - 8, and 

c = 3 X 10^° cm/sec have been used in relating t and T.] For = 1 hr., 

G = (2 X 10"^ V/m)^, values which we feel appropriate for the present model, 

T == 1 unit thus corresponds to approximately 18 hrs. Also, for these values the 

diffusion coefficient Dj^j^ = 1.5 x lO""^ L'’ Rg^/day. 

Note that the peak value of <Q > moves toward decreasing radius as time 
progresses. Asymptotic expressions for the position of maximum <CQ> , 
are derived in Appendix B. For L . = 8 and in terms of the parameter T, these 
expressions are 

L„.x = 8 [i-.42T + 0(t=)] T«1, (25a) 

" (ST)!/ “I T" T»1 . ^25b) 

For each value of T in Figure 1, the position of is accurately predicted by 
Equation 25b. 

Note further that the drop-off of <(Ql>with L is more rapid on the low L side 
of the maximiun than on the high L side. This feature is a consequence of the 
fact that the particle E xB drift velocity becomes progressively slower (owing 
to its 1/B dependence) as the dipole magnetic field source at L = 0 is approached: 
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particles on the leading inner edge of <( q > have the most difficulty executing 
radial drifts. 

On the other hand, our model is quite unrealistic at large L, where the 
combination of a rapidly decreasing magnetic field and a radially constant 
(though d, fp dependent) electric field yields unboimded electric drifts. As a 
result of such drifts there is a finite flux of particles to L ~ and the total 
number of M, J particles at all a is not conserved; i.e., from Equations 20 and 24 


dr 


^0 


<Q>da 




a<Q> 

dL 


03 


0 



I 

7777 



r>+ 1 

dx (i - x“ ) < 0 . (26) 

•'-1 


Aware of, but untroubled by, this fact we feel that proper modification of the 
model at large L would result in a particle conserving <Q]> essentially the same 
as Equation 24 at L->values of interest. 

Further feeling for the diffusion rate is gained by calculating the velocity 
with which moves inward. From Equations 25 we obtain 


= - 3.4 ^ [l+ 0 (T)] 

r , K, 

= - 4.8 X 10®T^Q[l + 0(t)J 


(T « 1) 


ft 


0.7x 10"<5 , 

« — zr—p, hrs 


\ c 




(27a) 
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max 

“"OF" 


2 flT 


.039 

1% 1 


[i + o(r^)] 


(T ^ ’> 1) 


r / «. 1 \ 

1 

f 0 . 7 x 10 ”®, \ 

[l 

y Hr 

hrsj 


(27b) 


In Equations 27 t and are again to be expressed in hours and U in(V/m)^. 

By way of comparison we have also solved the loss free diffusion equation 
for <Q > using the coefficient 


D, 


LL 


0.031 L, 


2 (±Y 

^ \k/ 


(earth radii 


3.6raO-Lj, [T 


7 T ^ (earth radii 


sec 


(28) 


calculated by Nakada and Mead (1965) for J ~ 0 particles diffusing under the 
effect of electric fields associated with sudden magnetic commencements. In 
this equation Ljj is the L value of the quiet time boundary at the sub-solar point 
in the Mead (1964) model of the magnetosphere (We take Lj, = 10.). The Nakada- 
Mead diffusion coefficient is used as a value rep'^eitentative of theories in which 
the driving electric fields arise from current variations on the magnetopause. 
For L^j = 10, the ratio of this coefficient to that 'a ssociated with a typical 
convection electric field is 2 x 10"® . 


Relating as given by luation 28 to via Equation 14, and using the 
definition L = - we find for this case 


D 


aa 


3.6 X 10 


-IS 


The diffusion equation in terms of normalized time r ' is thus 

c><Q> _ d / 1 a<Q>\ 


(29) 


dr' 


da 1^6 da j * 


(30) 
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For the same S -function initial condition and in the same notation used pre^dously, 
die solution to Equation 30 is 



In Figure 2 we display results of numerically evaluating Equation 31 with 
Lj = 8 as previously. *sQ]> is once more in arbitrary units with the same norma- 
lization as in Figure 1, so that direct comparison of the figure is possible. 

Note, however, that in this case T' -- 190 L.® 'R^ = .043 t (in hours), so 

that T' = 1 corresponds to approximately 23 hours of real time t . For values 
~ 1 hr., d - (2 X lO”"^ V/m)^, T' is thus roughly 4/3 times greater than T 
which parameterizes Figure 1. 

Qualitative features of Figure 2 are similar to those of Figure 1; the peal« 
in moves inward with time; <Q > drops off more sharply on its low L side 
than in its high L side; and there is a net flux of particles to L = <». The position 
of determined from an asymptotic analysis of Equation 31 similar to that 

carried out in Appendix B for Equation 24 is 


" 8{l- 1.6x T' +o[(T')^]} T' « 1 , 


(32a) 


max 


/ 


1 - ^ + 0 
T' 




T' » 1 . 


(32b) 
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The velocity of the <Q>peak in the two asymptotic regimes is 


dL, 


max 


■ar 


clT' 

- , 13 [l + 0(T' )] 


T' « 1 


R. 


5.6 K 10“^ [l + 0(t)] -p- t 23 hrs, 


(33a) 


dL 


max 


dt 


( . 3T' ) 








t » 23 hrs. 


(33b) 


As a final step we compare the velocities dl.^^jj^y'dt for the two mechanisms, 
calculating the ratio R of the rates predicted by Equations 27 (where the diffusing 
mechanism is the convection electric field) to the NaIcada~Mead rate Equations 33. 
ha the long and short time asymptotic limits R is respectively 


4.8 X 10 ^ T a 

R = — % 34 t « 18 hrs 

5.6 X 10 3 


(34a) 


R 


.039 ^ 1.5 

(r^a)'/n5/4 1-7 " 


t » 18 hrs 


(34fo) 


In computing R, values = 1 hr., fl = (2 x 10"'^ V/m)'^ have been used. 

Equations 34 express analytically a significant feature of Figures 1 and 2: 
at times when <Q> is sufficiently peaked (at an L value where diffusion is not 
inhibited by a strong magnetic field) , in the case of convection electric fields 
moves a good order of magnitude faster than in the situation where surface 
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currents are the agent. While Equation 34b does predict that the role of the 
current variation mechanism eventually exceeds that of convection electric fields, 
it is evident from the figures tliat this occurs only after has moved inward 
considerably and its velocity Of motion is measurably reduced from the early 
fast moving stage of its evolution. Compare also tlie height of the maxima in. 
Figures 1 and 2 for the same value of L; as time advances the maximum value 
of <!q >in Figure 2 is down from that in’ Figure 1 by over an order of magnitude. 

The features mentioned in the previous paragraph are the basis for our con- 
clusion that convection electric fields can move radiation trapped in the mag- 
netosphere radially inward at a rate at least an order of magnitude more rapidly 
than electric fields of magnetopause origin. 

CONCLUDING REMARKS 

While our results are encouraging in their prediction that convection electric 
fields can possibly be a strong diffusing force on trapped magneto spheric radia- 
tion, it would be hazardous to conclude that the mechanism presented here is the 
only agent. Sucii a conclusion is especially unwarranted in view of the small 
amount of evidence for either plasma flow or large scale electric fields in the 
magnetosphere. It may well be that a combination of effects are operative: 
diffusion by magnetic variations complementing the convection mechanism near 
t^e magnetopause and diffusion due to drift instabilities and even non-adiabatic 
processes playing a prominenc part at plasmapause distances and closer in. 
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In view of the scarcity of electric field measurements, an enlightening 
exioerimental exercise is to analyze statistically plasma flow and magnetic field 
in the magneto sheath. The autocorrelation time and mean square deviation for 
both these quantities is information from which corresponding properties of the 
magnetospheric electric field might be inferred. 

On the theoretical side, the addition of plausible loss mechanisms and so- 
lution of either the ensuing equilibrium or time dependent transport equations 
are needed before comparison with experimental flux'profiles becomes signi- 
ficant. Depending upon the L and <Q)‘ variation of t^e -diffusion coefficients, 
different solutions to the transport equation are expected. Comparison of such 
results with experiment may then further enhance the credibility of particular 
mechanisms both for inwardly diffusing particles and scattering them into the 
loss cone. 
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APPENDIX A 


DERIVATION OF THE a DIFFUSION EQUATION 


A derivation of the one dimensional diffusion Equation 5 is presented 
here. Our starting point is the two dimensional diffusion equation 

■3<Q>(a. A M, J, t) , a<Q> , / '5<0> 

5t ^ <“> ~W 


-1 

da 



~air] 




3<Q> \ 

5a / 





(A-1) 


derived by Birmingham, Northrop, and Falthammar (1967). Equation A-1 de- 
sci’ibes the evolution of the density (q) in an a, /S, M, J phase space, of particles 
acted upon by an ensemble of electric and magnetic forces each realization of 
which varies in time in an erratic fashion. The density (q) has been averaged 
over this ensemble of fields. Effects of the driving forces are included in the 
ensemble averaged velocities {a) and (yS) as well as in the diffusion coefficients 
^aa» ®/3yS’ S^^ce M and J are conserved during interaction with the 

fields, diffusion only occurs in the subspace of the Euler potentials a and 
For the dipole magnetic field we choose a = - /x sin^ i^r and /3 = 0. 
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To obtain the one dimensional a-equation, we first average (A~l) over 
at fixed a, M, and J. Denoting this ^average by a bar (~), we thus obtain 


d(Q)(a, M. J, t) 



(q>) 




+ 


d 




(A-2) 


Here use has been made of 1) the periodicity in /9 of all physical quantities and 
2) the canonical nature, 


da 



0 , 


(A-3) 


of a and /3 as guiding center coordinates (Northrop, 1963). 

Equation A-1 is correct through o(eS^) in the 0(e = m/e) adiabatic ex- 
pansion parameter and the 0(eS) (e « S « 1) smallness of the driving fields. 

We next assume that (q ) deviates from /3 -homogeneity only by an amount pro- 
portional to 8, i.e., AQ = (q) ” (Q ) = 0(8). This situation is satisfied, for 
example, if {q) (t = 0) is /3 -independent, since inhomogeneities in /5 then appear 
only as a consequence of the 0(eS) fields acting for a time (~e"^)of the order of 
a particle drift. As a result of this assumed near homogeneity in/S, Equa- 
tion A-2 simplifies to 
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To obtain the desired result, Equation 5, we must now show that 




(A-5) 


is of a negligible order of smallness. The canonical nature (Northrop, 1963) of 
the a, 0 coordinates is next used to show 



c dK(a, 0, J, M, t) 
e 


(A-6) 


For our model in which no magnetic variations occur, the Hamiltonian K is the 


total particle energy, kinetic plus potential. 

We are thus left with indicating the smallness of d /da A^a ) a(q)] . If in 

* 

Equation 2 for the electric potential, (a) = 0, we could straightforwardly argue 
that A {a), then at most o(e^ 8^ ) (no (a) arises from the dipole field alone) 
coupled with AQ ~ 0(8) leads to a negligibly small o(e^ 5^) term. The assump- 
tion that {a) = 0, however, means that the probability of a dawn-dusk directed 
equatorial electric field is the same as a dusk-dawn directed one. We feel that 
such equal probability is in contradiction with the extant circumstantial evidence 
and hence we here attempt to prove the smallness of d/da A{a)A(Q)j without 
the benefit of {a) = 0. 

Using Equations A-1 and A-4 it is readily shown that 


a(A(Q)) 

at 


+ CU. 


■9 (a<Q>) . _ \ ^<Q) 

— 5 ^ - aa 


0^e8^) , 


(A-7) 


where is the longitudinal drift rate in the absence of E. Equation A-7 has the 
solution 


a<q) 




t, M, J) 



- t), M, J, 



(a, M, J, f) , 


(A-8) 


i//(a,/3,M, J), the initial value of A(q), being heneeforth taken equal to zero. 
Employing the relationship (BNF, 1967) 



aA<v) 


(A-9) 


and the ejjplicit form Equation 2 for the potential, we find 
A (a) Aq] = 


{ a ) ^ fj?’ 
2 



dt' cos cOq (t' - t) 


a(Q> 

da 


(a, M, J, 



(A-10) 


(a) has been assumed to be independent of time. 

The presence of the factor cos ( t' - 1) oscillating over the drift period 

together with the slowly varying d {Q}/da render ^/^aTA^a) AQjnegligible as 


a driving term in Equation A-4. This may be seen formally by averaging 


Equation A-10 over a drift period ; 





Aq] 


<a)V^ ± fl. d ^ <Q > 

^ (27r)* do. d't 



(A-U) 


where o('^d/'^ci) correction terms in (A-11) are small in the ratio of the drift 
time 'Tjj (an time) to the diffusion time (an e“H“^ time). The loading 
order term in (A-11) is itself of o(rp/r^) in smallness compared with the dif- 
fusion term on the right hand side of Equation 5. Our conclusion is then that 

all terms represented by Equation A-11 are of a negligible order of 

) 

smallness and that (A-4) may be written 


^(q> J_ (fr- . 

dt da \^aa da j 


(A-12) 
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APPENDIX B 


THE ASYMPTOTIC MOTION OF L. 


We here calculate asymptotically as a function of time the L value at which 
(Q ), as given by Equation 24, maximizes. An analogous procedure has ^Iso 
been applied to Equation 31 but this calculation is not detailed here. 

For present considerations normalization of (Q ) is unimportant and Equa- 


tion 24 is written 


Cq> 


L3 rV4 


h ^ ± 

16 R/ r \Li‘» 


^ dx ( 1 “ exp 


SRjfL.^L^ 


(B-1) 


^max determined from the equation 




J_ 

4Rg'»r 


4Re'‘L.2l2t 


r 4 

I dx X (l - x^ exp — o~ 

^ 4 

I dx (l - x^ ) ^ exp T — ^ ^ — 

J_i ^ ' SR/L.^L^r 


= 0 . (B-2) 


We see no way to obtain analytically an exact solution to Equation B-2. However, 


progress is possible in the limits j3 - L? >> l, «: i when the 
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third term in (B-2) , 


ft* 


I = - 2(Sjj^\ln 


f dx^ (l exp /3x 

•'-1 


2ft {In G) 


(B-3) 


can be asymptotically evaluated. 

Consider first the short time r 0, /3 limit. In this case we introduce 
the change of variable x = 1 -y and write 


G = exp {ft) I* dy (2 “ 

Jo 


exp - fty 


(B-4) 


As the only contribution to the y -integral comes from the region about 
y = 0, where we expand 

(2-y)>/^ = (2)*'"' [l-| + 0(y2)] . 

Substituting this expansion into Equation B-4, we perform the integrations and 
obtain 

G = [t(t’ 2^) - W ^ °(^’^)] • (B-6) 
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Here r(a, x) is the incomplete gamma function (Gradshteyn and Ryzhilc, X965). 
For large x, y has the asymptotic form 


y(a, x) == P(a) + exp - x) (x -* oo) 


(B-7) 


where F is the ordinary gamma function. In the limit ,^5 G thus becomes 


. ( 2 )‘''‘>exp^ r / 5 \ 1 / 9 'i 

° ■ ^s74 " -giS 1 \?/ +0(/^ 2) 


(B-8) 


and 


I = - 2/3 


1 “ "J? + ■' o{6~^) 




(B-9) 


Substituting (B-9) into (B-2) and identifying /3, we finally write down the equation 


3 + 


4Rjfr\L^ L?L2 


c /8 L 2 L2 

+ -o + of— r 1=0. 


(B-10) 


The asymptotic solution to Equation B-10 is 






J 


(B-11) 


in the limit 1 ^ 1. 

In the opposite, long time, asymptotic limit, r-co, /3 -• 0, we approximate G by 


=-- £‘ dx ( 1 - x 2)V4 [l ^ 


(B-12) 
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Contributions arise only from eyen powers of x in the series. Performing the 
integrations, we find 


( 2 ) 




(B-13) 


so that in this limit 


I = - ^ + 0(/3‘>) 


(B-14) 


Lmax is hence determined from the equation 


3 + 


1 1 




4 RJ* r 


^ Rjf r2 L .^ 


+ 0 






= 0 , 


(B-15) 


whose solution yields 


A^ 


m a X 


EE(12r)>/‘* 


1 - 


112Re‘^L/t 


+ 0 
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_^Y 


(B-16) 
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FIGURE CAPTIONS 

<Q) , as given by Equation 24, for a ^-function at L. = 8 at T ?= 0, 
Diffusion is driven by convection electric fields, and for representa- 
tive magneto spheric conditions T = 1 unit corresponds to 18 hrs. 
of real time. 

{q ) from Equation 31 for an L. = 8, S -function initial condition. The 
diffusion coefficient is in this case taken from the work of Nakada 
and Mead (1965). T' = 1 unit corresponds to 23 hrs. of real time. 




Figure 2 


